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NOMENCLATURE 


A 



Fi 

k 

l 

M 

P 

Px 

Pr 

r 

r ef f 

r j 

T 


u,U 

uv 

V 

v j 

W 

w(C) 

Le 


P t pur 2 /ip 

afterbody (pressure) drag coefficient 

afterbody pressure coefficient, C p - ( p - p e )/ip e u e 2 

specific heat 

maximum afterbody diameter 
Xi/W 

enthalpy of species 
turbulent kinetic energy 

length scale for Prandtl mixing length model 

Mach number 

pressure 

pressure gradient in x direction 
turbulent Prandtl number 
radial distance from axis 

radial distance to effective inviscid boundary 

nozzle exit radius 

static temperature 

frictional velocity, u* * T w /pu 2 

axial velocity 

proportional to turbulent shear stress 

radial velocity 

injection velocity 

mixture molecular weight 

Cole’s universal wake function, Eq . (15) 

turbulent Lewis number 

net rate of production of the i th species 

vi 



X 


axial distance 


core length 

mole fraction of i^ species 
distance normal to inviscid plume interface 
normal distance to point where u ■ £(112 + 1 ^) 
ratio of specific heats 

boundary layer thickness; also characteristic mixing layer thickness 
boundary layer displacement thickness 
turbulent dissipation rate 
turbulent viscosity 
density 

shear layer spreading parameter, defined by Eq. (44) 
spreading parameter for U 2 /u-^ * 0 
Prandtl number for turbulent kinetic energy 
Prandtl number for turbulent dissipation 
wall shear stress 
stream function 
Subscripts 
c centerline 

e,E external flow 

j , J exhaust plume (jet) flow 

1 inner mixing layer boundary 

2 outer mixing layer boundary 

t stagnation 

I inner stream 

i refers to i th species 
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1. SUMMARY 


The development of a computational model (BOAT) for calculating nearfield 
jet entrainment, and its incorporation in an existing methodology for the 
prediction of nozzle boattail pressures, is discussed. BOAT accounts for the 
detailed turbulence and thermochemical processes occurring in the mixing 
layer formed between a jet exhaust and surrounding external stream while 
interfacing with the inviscid exhaust and external flowfield regions in an 
overlaid, interactive manner. The ability of the BOAT model to analyze simple 
free shear flows is assessed by detailed comparisons with fundamental labora- 
tory data. The overlaid procedure for incorporating variable pressures into 
BOAT and the entrainment correction employed to yield an "effective" plume 
boundary for the inviscid external flow are demonstrated. This is accom- 
plished via application of BOAT in conjunction with the codes comprising the 
NASA/LRC patched viscous/inviscid methodology for determining nozzle boattail 
drag for subsonic/transonic external flows. An assessment of the overall 
approach is provided via comparisons between these predictions and data for 
underexpanded laboratory cold air jets. A manual describing the operation of 
the BOAT code comprises a supplement to this report. 

2 . INTRODUCTION 

The accurate prediction of nozzle boattail drag requires a detailed 
description of the coupled viscous/inviscid flow processes occurring along 
the nozzle afterbody and in the nearfield mixing layer growing along the 
plume interface which separates the nozzle exhaust and external air streams 
(Fig. 2.1). While analyses based upon solving the turbulent Navier-Stokes 
equations can provide a description of this complex flowfield, the widely 
disparate length scales and flow characteristics in the various regions in- 
volved lead to rather prohibitive computer time requirements in achieving 
results of adequate resolution. A more efficient procedure is that provided 
by a "patched" viscous/inviscid methodology where each region is separately 
analyzed by computational procedures specifically catered to the flow pro 
cesses and length scales occurring within that region. The overall flowfield 
solution is arrived at by patching these regional solutions together in an 
iterative manner. Such a patched methodology has been implemented at NASA/ 
LRC*’ 2 ’ 3 for subsonic/transonic external flows. For nonseparated afterbody 
flows, the LRC system presently includes the relaxation procedure of South 
and Jameson 4 for analyzing the inviscid subsonic/transonic external flow, 
the supersonic exhaust plume model of Salas , and an extended version o t e 
Reshotko-Tucker 6 turbulent boundary layer integral method. 




Figure 2.1 Schematic of afterbody/jet exhaust flowfield. 

In past applications of the patched LRC system to af terbody/exhaust flow 
configurations , the inviscid plume interface was treated as a solid surface 
in determining afterbody pressure distributions. This assumption yielded re- 
sults which substantially underpredicted the afterbody drag, since the addi- 
tional contribution of jet entrainment due to mixing processes along the 
interface was not included. The physical influence of jet entrainment is 
manifested in an overall reduction of the blockage effect of the inviscid 
plume. The incorporation of entrainment effects into the LRC system by a 
displacement thickness type correction to the inviscid plume boundary will be 
shown in this report to yield substantially improved drag predictions. 

The BOAT code, developed to calculate these entrainment effects, solves 
the axisymmetric jet mixing equations for a reacting gas mixture by a mixed 
implicit/explicit finite-difference procedure in transformed (x,i(i) coordinates. 
In developing BOAT, the "best" features of several current mixing/afterburning 
models were combined to achieve rather unique computational capabilities. In 
this respect, BOAT employs: (1) the same generalized treatment of thermo- 

chemical processes and implicit/explicit streamline integration procedure used 
in the Low Altitude Plume Program (LAPP) 7 (the current JANNAF standardized 
model for analyzing mixing/afterburning processes in low altitude rocket 
exhaust plumes); (2) the same type of shear layer discretization, grid distri- 
bution, and growth rules used in the GENMIX code 8 ; and (3) the same type of 
overlaid procedure for describing variable edge conditions and both normal and 
streamwise pressure gradients as used in the GASL "patched" system 8 for the 
detailed analysis of inviscid/shock and mixing/afterburning processes in rocket 
exhaust plumes. Previous entrainment models* 8 ’** have not accounted for 
pressure gradients and thermochemical effects. The present work will show 
that the neglect of pressure gradients can cause serious errors in the pre- 
diction of nearfield entrainment. While no detailed calculations have as yet 
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been made for hot, afterburning plumes, it is apparent that real gas behavior 
and chemistry will affect nearfield entrainment. 

Of critical importance in the overall modeling of nearfield jet entrain- 
ment is the selection of an appropriate turbulence model. Several turbulence 
models of varying complexity are presently incorporated in BOAT. For the low 
speed, nonreactive cases thus far considered, an extended version of the 
Prandtl mixing length (ML) model and/or kc2 two-equation turbulence model 
appear to yield results of sufficient quality. A description of these models 
and their assessment via comparisons with pertinent free shear layer experi- 
mental data will be discussed. 


Applications of BOAT to the analysis of laboratory cold air jets have 
been performed and sensitivities to various parameters explored. Detailed 
experimental measurements in the underexpanded jet nearfield region are not 
presently available to adequately assess the individual modeling procedures 
employed. The validity of the overall approach must therefore be establishe 
by comparisons between predicted boattail pressure distributions (via appli 
cation of BOAT in conjunction with the other components of the patched LKL 
system) and measured boattail pressures for the cold air jets . These com- 
parisons are presented herein and the sensitivity of predicted boattail pres- 
sures to various model parameters is discussed in some detail. 

The overall methodology is discussed in Section 3, while the computa- 
tional procedures incorporated in the BOAT code are detailed in Section 4 of 
this report. Section 5 describes the validation of the code via comparisons 
with data used in the 1972 NASA Langley Free Turbulent Shear Flows Confer- 
ence 13 . Section 6 describes its application to laboratory cold air jets, 
which simulate conditions for actual aircraft nozzles. The supplement to 
this report contains a detailed description of the operational procedures 
employed within BOAT, including a complete listing and comprehensive users 
manual. The displacement thickness-type correction introduced in this stu y 
and the methodology for incorporating BOAT in the overall Patched system were 
lointly formulated by the authors and Dr. R.G. Wilmoth of LRC, who addition- 
ally performed the systematic calculations for the cold air jet cases re- 
ported in Section 6. 


3 . METHODOLOGY 


3.1 Background 

In developing the BOAT code, a review of several widely used, well tested 
models was undertaken with the objective of assessing their best features and 
incorporating these features into BOAT. The models reviewed included the LAPP 
code 7 , the GENMIX code 8 , and the GASL plume model . 

The popularity of the LAPP code 7 is largely attributable to its ease of 
usage and overall reliability. LAPP calculates the constant pressure mixing 
between concentric, chemically reacting streams. The governing parabo ic 
mixing layer equations are solved in Von Mises (x,*) coordinates using a mixed 
implicit/explicit computational procedure. In particular, the species 
tinuity equation is solved in a coupled fashion using an implicit, linearized 
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procedure which permits substantially larger integration step sizes in chem- 
ical near-equilibrium situations than explicit procedures. The computational 
scheme requires the definition of grid points from the jet axis to the upper 
mixing layer boundary, thus precluding a refined definition of nearfield 
shear layers as well as the treatment of initial boundary layers. Mixing 
growth occurs via the addition of grid points at the upper boundary 
which is numerically cumbersome and can be rather inaccurate. Arbitrary 
chemical reaction mechanisms can be treated via the input of appropriate 
forward rate constants and thermodynamic data. 

The BOAT code closely follows the overall methodology employed in LAPP. 
Since only the nearfield shape of the "effective" plume boundary (i.e., its 
shape in the first Inviscid cell or wave length) significantly influences 
nozzle boat tail pressures t the principal modifications in this methodology 
have been geared to improving the computational accuracy in the nearfield 
shear layer. A rather efficient approach for treating nearfield shear layers 
was introduced by Patankar and Spalding in the GENMIX code 8 . Here, the com- 
putational domain spans the shear layer which is solved in a mapped x,oj 
coordinate system*. The growth of the shear layer is determined by specified 
rules related to property variations at the edges. This approach has been 
adapted in BOAT. 

The transformation to x,u) coordinates in GENMIX numerically confines 
the shear layer between oj = 0 and 1, thus permitting its calculation to be 
performed with a fixed number of grid Intervals. The equations In x,w 
coordinates are lengthier than those In the x,\p system and the integration 
process does not follow streamlines (i.e«, lines of constant oj are not 
streamlines). The desirability of working in an x,ty system Is quite evi- 
dent. However, previous models such as LAPP have not implemented this 
system In an effective manner in nearfield situations. Patankar and Spalding 8 
discount the use of an x,^ system entirely stating that "advantageous 
though (its use) may be in other respects, the grid nodes are still ineffi- 
ciently distributed . " In the BOAT code, however, a new computational pro- 
cedure is employed which uses the x,ip system for integrating the flowfield 
equations, yet efficiently distributes the grid nodes In a manner identical 
to the mapped x,w system in GENMIX. ~ " 

The requirement for employing a discretized shear layer approach In 
studies keyed to the accurate prediction of nearfield entrainment is quite 
evident from the results depicted In Figure 3.1. Here, BOAT with a fixed 
11-point grid array is shown to predict the correct linear relation between 
mass entrainment and axial distance while LAPP, initialized with 21 points 
across the jet, requires more than two exit radii before the proper entrain- 
ment rate is asymptotically reached. In terms of profile shape, the LAPP 
prediction was poor at the two— radii location, showing a lack of convergence 
to the correct solution, while the BOAT calculation Is in excellent agreement 
with experimental data, as will be detailed in Section 5. 


In GENMIX, the mapping is given by u = & - *l ( x J / 0»2 (x) _ *1 (xl] where ^ 
and if>2 are th e streamfunction values at the upper (i^ 2 ) and lower (i|n) shear 
layer boundaries. 
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For the incorporation of inviscid structure into the mixing/afterburning 
analysis, the '’overlaid” approach introduced in the GASL plume model 9 has 
been adapted. The overlaid concept is a direct extension of classical 
boundary layer methodology to the analysis of nearfield shear layers. In 
boundary layer theory, the inviscid flow pattern is first calculated followed 
by a boundary layer calculation with edge conditions and pressure gradients 
set by the inviscid flow pattern. In the direct extension of this approach 
to nearfield shear layers, the inviscid exhaust plume and external flow 
patterns are first determined and the nearfield shear layer calculation is 
then initiated along the inviscid plume interface (see Fig. 3.2). Local edge 
conditions and both normal and streamwise pressure gradients are set in 
accordance with the calculated inviscid flow pattern and the rate of growth 
of the shear layer. The validity of this approach in application of the GASL 
model has been established by comparisons with various sets of experimental 
data. Selected comparisons are provided in the survey paper by Dash and 
Pergament 14 . 

X SUBSONIC /TRANSONIC EXTERNAL FLOW MAP 
° SUPERSONIC EXHAUST PLUME FLOW MAP 
• MIXING LAYER GRID 



Figure 3.2 Schematic of mixing layer overlaid on 
jet exhaust and external flow maps. 

3.2 Governing Equations 

BOAT solves the parabolic jet mixing equations in transformed (x,^) 
coordinates. The resultant system of equations is listed below and includes 
equations for the turbulent kinetic energy, k , and dissipation, z , as 
required in the two-equation turbulence model option. 
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9k 
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Turbulent Di ssipation 


li-il* (tH) + i( c i p - c 2 t ) 


where 


A = y 


phi; 

t 4 , 


Au / 3 uN 

4> V3 4>J 

and the transformation from (x,r) to (x,i|>) coordinates is given by 


(3) 


(4) 


(5) 


^ fr “ pUT 

31 ( 6 ) 
* H * - pvr 

The above equations adequately describe mixing/afterburning processes in 
the nearfield shear layer for mildly underexpanded plumes, as well as in the 
farfield mixing layer. For plumes with significant underexpansion, nearfield 
mixing processes are more appropriately described in plume-oriented boundary 
layer coordinates measured along and normal to the inviscid plume interface. 

A discussion of the equations and methodology entailed in this system is pro- 
vided in Reference 14. 
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3.3 Effective Plume Boundary Concept 

The mechanism for modifying the inviscid plume geometry to account for 
entrainment is a direct extension of the weak interaction approach in stan- 
dard boundary layer (BL) theory. In BL theory, the solid wall shape must 
be modified so that the "new" effective wall shape Induces streamline de- 
flections in the inviscid solution compatible with the viscous flow pattern. 
A standard modification involves applying an Injection velocity boundary 
c - d , ition - ’ a l° n S the solid wall shape, or adding a displacement 

thickness, 6 (x) , to the solid wall shape (Fig. 3.3). For two-dimensional 
flows and axisymmetric flow situations where the BL thickness is small 
relative to the transverse body dimension, the normal velocity at the BL 
edge, v e (x) , may simply be applied along the solid surface. For "thick" 
axisymmetric boundary layers, the v(x) distribution applied along the 
surface must be modified from the edge distribution, v e (x) , to satisfy 
the continuity equation. The displacement thickness variation, 6*(x) 

Is readily determined from the normal velocity distribution, v(x) , by 
solution of the differential equation 

d*5* _ v 

dx u (?) 



Figure 3.3 Effective boundary concept for a wall boundary layer. 

In extending this concept to shear layers and the resulting modification 
of the inviscid plume geometry, we simply seek an appropriate boundary con- 
dition for the inviscid external flow calculation that will reproduce the 
normal velocity distribution, v e (x) , at the outer edge of the shear layer 
(Fig. 3.4). If the entrainment induced flow deflections are small in compar- 
ison to the streamline deflections induced by the inviscid plume geometry 
(i.e., the "blockage effect"), the analogy with the BL problem is quite 
s ^ ra ^8^t f orvard . Referring to Figure 3.5, the inviscid plume geometry can be 
properly modified by prescribing the inflow velocity, v e (x) , along the 
inviscid plume interface in the repetition of the external flow calculation. 
Equivalently, the displacement thickness, 6*(x) , can be added to the 
inviscid interface yielding a new "effective" plume geometry. The specific 
choice between these two methods is somewhat arbitrary. The effective geom- 
etry concept provides a smoothing effect via the integration of the v e (x) 
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distribution and was chosen for use in the subsequent calculation primarily 
because the external flow code 1 * presently utilizes a geometric type boundary 
condition. 


U 


e 


u e 



Figure 3.4 Effective boundary concept for 
an axi symmetric shear layer. 



Figure 3.5 Addition of displacement thickness to inviscid plume 
interface yielding effective plume boundary. 

The validity of the effective boundary concept can be demonstrated by a 
simple example, illustrated in Figure 3.6. Here, the plume is perfectly ex- 
panded and a constant pressure mixing calculation is performed in the absence 
of initial boundary layers. The outer edge normal velocity variation, v e (x) , 
predicted by a BOAT calculation, was then integrated to determine 6 (x) [via 
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Eq « (7)] yielding the illustrated effective plume geometry (obtained by adding 
6 (x) to the plume interface r. «= 1.0 ) . Performance of the inviscid ex- 
ternal flow calculation over this geometry yielded a normal velocity distri- 
bution along the upper shear layer boundary which compared favorably with 
BOAT s predictions for v g (x) , Indicating a matching of the streamline slopes 
along this boundary. An accurate representation of the entrainment effect on 
the upstream inviscid flow is therefore expected. The slight disagreement at 
the initial shear layer region Is produced by the numerical smoothing (In the 
Inviscid code) of the effective shoulder In the boundary (I.e., the flow is 
forced to turn smoothly Instead of the sudden turn predicted by BOAT) . The 
slightly increasing difference between the Inviscid and BOAT-calculat ed normal 
velocities, with increasing x , is attributed to prematurely terminating the 
effective body at x/rj s 3 and employing a cylindrical boundary downstream. 




Figure 3.6 Comparison of effective entrainment velocity calculated 
with BOAT and with inviscid external flow code using 
effective boundary concept , u e ^ u j “ 0.2 , M e ■ 0.4 . 



4. COMPUTATIONAL PROCEDURES 


4.1 Integration Scheme 

In BOAT, a parabolic system of equations [Eqs. (l)-(5)] are integrated 
in the axial direction from specified initial conditions u(x-p\|0 , T(x-p40 * 

Fi(x lf *) , k(xi,*) , and c(x lt *) where x 1 is the initial axial station 

and w ^ ere ^i( x i) and ^2^ x l^ are t ^ e u PP er and lower mixing 

layer boundaries. A fixed number of grid points, (N) , are equally spread 
across the shear layer in intervals of A^> = (^2 “ ^i) / • 

Referring to Figure 4.1, the following steps are performed in integrating 
the flowfield equations from x^ to X 2 ( = + Ax) : 

o The rate of mass entrainment is determined at the upper and lower 
boundaries at x^ (see Section 4.6) and used to evaluate and 

^2 at x 2 in an explicit fashion [i.e., ^i f 2^ x 2) “ ^x^ 2 ^ x l^ x ^’ 

o The dependent variables at the initial station x x are redistributed 
over the extended interval * 1 (x 2 ) * * * ^ 2 ( x 2 > equally spaced 

Intervals of Aip = [^ 2 ( x 2^ ” ^1 ^ x 2^ ^ ^ ^ • 

o The equations are integrated along the streamlines *Kl) * ( x 2 ) 

+ (I - 1)/(N - 1) * txp 2 Cx 2 ) - 4^1 Cx 2 > ] for 1=2 to N-l , with edge 
conditions at 1=1 and I = N specified, consistent with the 
inviscid flow pattern. 

The sequential process above Is then applied to advance the equations another 
integration step, etc., with the first step above performed at the point in 
the program where the calculated variables at x 2 are reset into the initial 
profile array. 


4.2 Overlaid Procedure 

In analyzing und er expanded plumes, the detailed inviscid exhaust /ext ernal 
flow solution is used to specify the variable edge conditions and pressure 
gradients for the subsequent shear layer analysis. The inviscid solutions are 
supplied to BOAT via arrays at arbitrary axial stations for both the exhaust 
and external flows. BOAT processes these data yielding orderly vector arrays 
in mapped stream-function (x,<u) coordinates, which eliminates the need for 
inviscid map searching procedures. Referring to Figure 3.2, the procedure is 
as follows: 

o At axial stations X K , the user supplies arrays of rj , Up , 

P x > an d Tp . For the jet exhaust flow, these arrays extend from 
the axis to the plume interface while for the external flow, they 
extend from the interface to an arbitrary location. The axial 
stations for the jet exhaust and external flow need not coincide. 

o The program calculates the value of at each point I via 

Integration of Eq. (6) . For the jet exhaust data we have 
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while for the external flow we obtain 


♦i = 



(8b) 


o 


where ’f'jet is a constant and rj e t is the local height of the 
plume interface at the axial station of interest. 

The value of nondimensional stream function u> is determined at 
each point I . For the jet exhaust, cu is given by 

“I " Wt 


while for the external flow, by 

u>I = [*i-*jet 1 /I *ext (x> -*jetl 


(9b) 


o 


where <^ ext (x) is the variable value of 1(1 at the external flow 
upper data boundary. 

Vector arrays, equally spaced in us (0 - id i 1) , are created for the 
jet exhaust and external flows of the form VJET M>L>K and VE ^M,L,K ’ 
M designates the variable type (1 = r , 2 = P, 3 = T , and 4 = • ) , 

L designates the value of us [ us - (L - 1) / (Ljnax “ 1) where Lmax s 

a fixed integer constant specified by the user for defining the grid 
definition in the mapped arrays] , and I designates the axial 
station . 


Since all jet or external profiles have the same number of mapped data 
points, L, a v , radial map search procedures are also eliminated. Thus, 
properties It the lower shear layer boundary * (x) are determined as 
follows : 


o 


w * ^i(x)/i|>jet 
such that oj(L) 


is determined yielding the index 
- 0) - d (L + 1) . 


L 


w * (^nax + ^ 



o Properties ftf,K are determined at the stations K and R+l (where 

V— £ V S Y _ \ V^T f-Vs r-s 1 -J — _ J "1 _ » r 


x ^ X K+1^ th e li near interpolation 


M,K 
K + 1 


VECT M f L,K + ^KLH iK 


M,L,K 

K+l 


(10a) 


where *= a) * (Lmax * 1) * (D-l) and VECT represents either VJET 
or VEXT . 


o The desired properties f M are determined by the interpolation 

f H' f H,K + ra ( f «, K+ l-V) 


(10b) 


where RX - (x - x K ) /(x K+ 1 - x K ) 


For the external flow, the same approach is followed with w in the first 
step above evaluated by the relation u = [* 2 (x) - * iet ] / [+ ext (x) - et ] . It 
has been found most expeditious to evaluate the pressure gradients across the 
shear layer directly from the mapped arrays (i.e., p x (x,^ L ) = [VECT (2, L, K+l) 

- VECT (2, L, K)]/(x r+1 - x K ) where x R <x<x K+1 and to, = ♦t/iI'w } for each 
individual streamline, t|» l . . ^ Jei: 

4.3 Finite Difference Formulation 

A. 3.1 Eq uations . - The same implicit/explicit formulation employed in the 
LAPP code has been retained in BOAT. With the notation depicted in Figure 4.2 


m + 1 


m- I 



A \j/ 


A \js 


b Ax b 

n n + l 


Figure 4.2 Finite difference index notation. 
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the resulting difference equations are written: 
Species Diffusion 

( F i)n+l,m = ( F i)n,tn + 
and on the axis (m «= 1) by 


Le Hpu rj, ( \ 

Pr ^ \ i/'l' 


(*i)n+l,m tx 


(pu) 


tl,TD 


n,m 


( F i)» + l,l • ( F l)n,l + 2ix ( 1 fr\ >1 ** 


( F l)< 


i)n 4- l,m 


Ax 


n,l 


(pu) 


n,l 


Axial Momentum 


U n + 1 ,m U n,m ^ l ^ 

\ / n,m 


Ax 9__ ( 


ypur* 


^ P x^n 


Ax 


and on the axis (m - 1) by 


u n + 1,1 “ u n,l + 2Ax(y) n,l ^ W 


( Px ) n+ 2 >: 


Ax 


(pu) 


n,l 


n,l 


Energy 

T n + l,m ' T n,m + iFc~ 3^ iPr ^ *\pj 


Ax 3 _ (1 e hPur 2 T 


m p 


, ( Px ) n+ 2» 


n ,m 

Ax 
m 


n,m 


Ax 


( pc p)n,m ( puc p)n,m i 


■Z( h i*i)n,i 


v 'n.m m 


and on the axis (m = 1) by 


3^Wn,m + (pr) Wn,ffiX( C Pi) ( 

n ’ m i \ / n,m 




n,m 


(11a) 


(lib) 


(12a) 


(12b) 


(13a) 
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T , , , * T 
n+ 1,1 


«.i + 2 (£) n>1 h ( T *)„,i 

( P x)t,+ 1,1 4j< J 


2( h i 4 i)n,l 


( puc p) n ,l 

where first derivatives f ^ are given by 

( f \ _ ^n,m + 1 ^n,m - 1 

. Vn,m 2 A if) 

and second derivatives of the form ^af^ are given by 
|y (afA = a ( f n ? ni+l ~ f n,m) 

W V *Jn,m n>m+ l 

(f - f ) 
\ n ,m n,m - 1/ 

1 9 


- a 


n * m * 2 A 


where a . «= 

n,m±- 2 


The pressure gradient (p j 1 is given by 

' X/ r, J 


n + y , m 


(13b) 


(14a) 


(14b) 



P 


n + 1 ,m 


Ax 


P 

n,m 


(15) 


The above equations are of first order accuracy in Ax and second order 
accuracy in AiJ> . Analogous expressions result for the turbulence equations 
for k and e . 

4.3.2 Linearization of chemistry terms .— The species mole fractions at 
station n + l,m are determined from the species conservation equations by 
linearizing the chemistry terms [i.e., (^i) n + 1 m 1 and inverting the 

resulting matrix. The linearizations involving species and Fj (for a 

two-body reaction) or , Fj , and F^ (for a three-body reaction) at 
station n + 1 (all variables are known at station n ) are given by 
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/ \ / \ 

)„♦! ■ ( F i F j) n + F :n[( F i)n + l - ( F ijn_ 

+ F . 
in 

( F j)„ + 1 - \j/n 


■ - (Vj) £ + ( F j)n( F l)n + l + ( F l),( F j).+ l 

( F i F j F k)n + i ‘ - 2 ( F l F j F k)„ + (Vk)i( F i)„ + l 

+ ( F l F k) n ( F j) n + 1 + ( F i F j)„ ( F *)„ + l 


The terms underscored by a single line contribute to the elements of the 
coefficient matrix, while the terms underscored by a double line contribute 
to the known column matrix on the right-hand side of the matrix equation for 
the linearized system. Thus, the matrix equation takes the form (for N 
species) : 


a ll a 12 a 13 

3 IN 


_F l“ 


" Q l" 

a 21 



F 2 


Q 2 

®N1 

a NN 


f n 




(17) 


4.3.3 Allowable step size .- The maximum allowable marching step must 
satisfy explicit stability considerations for a parabolic system of equations 
as well as those of overall accuracy. As such, it is taken to be the 
minimum value of 

to.f- a* 2 (18.) 

\ 2Leppur z / 

' 1 n,m 


and 


Ax 


y 


n ,m 


y 


n,m 


1 


(18b) 


as evaluated for all grid points m at station n . On the axis (m 1) , 
Eq. (18a) is replaced by 


Ax 



(19) 


17 



The chemistry terms impose no stability limitation on the step size due 
to their implicit formulation. To insure an accurate description of chemical 
processes when the chemistry is "fast" (i.e., when one or more reactions are 
near equilibrium), a maximum allowable temperature change, AT max , is per- 
mitted along a streamline in an integration step. The value of ^T max is 
input to the code. Should the temperature change exceed this value, the 
step size is repeatedly halved until the temperature change is less than 
^max > or the step size becomes less than a user input minimal step, in 
which case the program terminates. Should the user input too large a value 
ATjaax for the case considered, the computed species mole fractions can 
become negative in "fast" chemistry regions. This is also handled by step 
size halving until the predicted mole fractions are all positive. The choice 
A T max = 5°K should suffice in the most stringent situations, but may be 
overly limiting in general (i.e., with a coarse grid, such changes could occur 
in a given step due solely to diffusive processes). In nonreactive situations, 
or situations where rapid burning is not anticipated, a value of AT max » 25°K 
may be more appropriate. 


4.4 Chemical Reaction Rate Equations 


Ten possible reaction types are included in the program: 


Reaction Type 


(1) 

A + 

(2) 

A + 

(3) 

A + 

(4) 

A + 

(5) 

A + 

(6) 

A + 

(7) 

A + 

(8) 

A + 

(9) 

A + 

(10) 

A + 


B Z C + D 

B + M Z C + M 

B Z C + D + E 

B Z C 

H Z C + D + M 

B + C + D 

B + M -► C + M 

B + C + D + E 

B ■+ C 

M + C + D + M 


Reaction types (6) - (10) correspond to reaction types (l)-(5), but proceed in 
the forward direction only. In reactions (2), (5), (7), and (10), M is an 
arbitrary third body. In this program, all species are assumed to have equal 
third body efficiencies; thus, in evaluating w^j) , = (W)”^ . The 

formulation for the net rates of production for each reaction type is written 
below. 


( 1 ) 

( 2 ) 


*0) - k f p 2 F F - ^- F C F P 

K f p Vb k 


(j) vVb V 2f c 


w 


K WRT 
P 
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, k p J F F F RT 

(3) w - k f p F F v 


.(j) 2 _ k f pF C 

w *= k f p F a F b K rt 

P 


• a) V^!a k f p3r c F p RT 

w = W K W 

P 


(6) v (J) - k f P 2 F/B 


(7) v 


0).vVb 


(8) w (j) = k f p 2 F A F B 


(9) w (j) = k f p 2 F A F B 


k f p 2 F 

(10) w (J) = ' 


forward rate coefficient, k f , is expressed in the form, 


k f *= AT _N exp (B/RT) 


and K p is determined from 


In K p = - AG/RT 
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The rate coefficients are divided into eight types: 


Rate Coefficient Type 

(1) * A 

(2) k f = AT -1 

(3) k f « AT -2 

_ 1 

(4) k f *= AT~ 2 

(3) k f = A exp (B/RT) 

(6) k f ■= AT -1 exp (B/RT) 

_ 3 

(7) k f = AT 2 

(8) k f « AT _N exp (B/RT) 


4.5 Turbulence Model 

4*5*1 Mixing l eng th model .- In the extended ML model employed in BOAT, 
the turbulent viscosity, p t , at all interior points is given by the 
expression 


y t 


p £ 2 


5u 

3y 


and on the axis by 


^t 


= pi 3 


9 2 u 

9y 2 


(22a) 


(22b) 


In application of this model, the mixing length l is linearly related to 
the mixing layer thickness 6 , i.e., 1/6 * constant . Several definitions 

of 6 are illustrated in Figure 4.3. We have attempted to achieve some 
degree of generality by the use of just two scaling constants (a nearfield, 
two-dimensional. constant and a farfield, axisymmetr ic constant). No further 
empirical relations are employed. The scaling constants, 1/6 , were 
determined by carefully matching the fundamental experimental data to be 
discussed in Section 5. 
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In application of this model in x,^ coordinates, at an interior point 
n,m we employ 



while at the axis, the relation 



Is employed. 


(23a) 


(23b) 


4*5*2 Two-equation model* , — The ke2 model, developed by the group at the 
Imperial College 1 5 , has been incorporated in the BOAT code. Here, the turbulent 
viscosity is given by 



(24) 


where Cy and the constants , C 2 , a k 5 a ^d a £ have been extracted 

directly from Reference 15. As in the ML model, the constants differ for 
nearfield (2D) and farfield (axisymmetric) situations. These constants are 
summarized below: 


C p - 0.09 g (P/e) - 0.0534 f 

c i - I-* 

C 2 = 1.94 - .1336 f 

o k = 1.0 

a e * 1 - 3 


h/2 

du c 

_ ^£~ 

_ u c 

dx 

dx 


where P/e ( P is the turbulence production rate) is the shear stress weigh ted 
average across the mixing region. The functional dependence of g on p/ E 
and of f on the axial velocity centerline decay follow Reference 15 directly. 
In this regard it is important to note that the axisymmetric correction term, 

* ' — 

The incorporation of the ke2 turbulence model into the BOAT code and its 
subsequent validation by comparison with experimental data (see Section 5) 
was supported by the Naval Weapons Center under Contract No. N00123-78-C-0010 . 
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f , was obtained by matching data on constant pressure, nonreacting coaxial 
jets. For generalized plume studies where the centerline velocity addition- 
ally changes due to both pressure gradients and chemistry effects, ^it may be 
necessary to isolate the diffusive decay effect by use of a "dummy" inert 
species. Several new approaches 16 ’ 17 (which do not depend on the centerline 
decay) for extending kc2 models to axisymmetric flows are presently being 
investigated and should lead to a greater degree of generality in application 
of this model to complex flow situations. 

A. 6 Shear Layer Growth 

The smooth growth of the shear layer is provided for by the use of rules 
dependent upon the profile variations near the edges. A certain amount of 
trial and error has been exercised in arriving at rules which are generally 
satisfactory. Following the approach of Patankar and Spalding , we start with 
the limiting expression 


m ■ fcim — (uty u y) v / u y 

Uy + 0 7 7 


(25) 


for the local rate of mass entrained at the shear layer edges and 

With the velocity variation near the edges approximated by the expression 


u(y) - u 1>2 + c(y- yi y 2> 

and using the mixing length expression for p t we obtain 
m n>1 = Ap n>1 i 2 (u n> 2 - u n) i)/(y n ,2 - yn,l) 2 

at the lower boundary, m = 1 , and 

»n,N = 4p n) N* 2 < u n,N" u n,N- l )/( yn,N-Fn,N- 1>‘ 


(26) 


(27a) 


(27b) 


at the upper boundary, m - K . These relations are applied in an explicit 
fashion in predicting the values of the stream function at x + Ax 

Thus , 


^l(x+Ax) ^l(x) + \ 

n , r 

and 

^2 (x + Ax) * ^2 (x) + ("^") njN AX 


(28a) 


(28b) 
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For the ke2 model, the shear layer growth at the upper boundary, , is 

estimated by the expression ^ 


+ Ax) = ( x ) + 3r n 


2 v t r 


/ r N r N - 2 \ 
1 \ r N ~ r N - 1 / 


(29) 


with a similar expression employed in estimating the growth rate at the lower 
Ol^) boundary. 


4.7 Initialization Procedures 

Three options are provided in BOAT for initiating the mixing layer calcu- 
lation; namely, a generalized option where the user specifies initial profiles, 
a shear layer option, and a boundary layer option. 

4*7. i Generalized option .- In this mode, schematized in Figure 4. A, the 
user inputs values of the dependent variables u and T , and the mole 
fractions* X ± , as a function of radial distance from the axis at an arbi- 
trary x initial station. The point 1=1 Is either the lower edge of the 
shear layer, ^i( x ) , or the axis, while the point I = N is the upper mixing 
layer boundary. In using this option in conjunction with the overlaid proce- 
dure, the edge conditions (at 1=1 and N ) must be consistent with the 
inviscid flow solutions as defined by the exhaust and external flow maps. 

After reading in the data indicated in Figure A. A, the mass flow variable 
is determined at all points I and an equally spaced profile in A^ is gen- 
erated, where 


A<f> = (4 > n - ♦j) / (N - 1) (30) 

Thus, the N specified data points, Input In arbitrary intervals in the phys- 
ical plane, are recast into an equally spaced array of N points in the 
transformed plane. The same number of grid points, N , is retained through- 
out the overall calculation. 

The generalized option should always be implemented when the initial pro- 
files are known. This, however, is not usually the case, and the options 
presented below can be utilized to yield code generated initial profiles. The 
code also has a restart option which allows for marching to a specified axial 
station and restarting from that point with no additional input requirements, 

A. 7. 2 Shear lay er op ti on .- In this mode, schematized in Figure A. 5, a 
fully developed turbulent shear layer is assumed to exist at the starting 
station. For use in an overlaid approach, the user need specify only the com- 
position of the jet and external streams, and the starting location x . Local 

* 

It should be noted that while the user inputs the composition in terms of 
, the integration is performed in terms of the dependent variable 
where F^ = X^/W and W is the local mixture molecular weight. 
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edge properties will be determined from the inviscid data maps. In the con- 
stant pressure mode, the velocity and temperature of the jet and external 
streams will additionally be required. 


For most situations of interest, the initial shear layer width can be 
reasonably estimated by the incompressible relation 18 



(31) 


where 1 and 2 designate the lower and upper shear layer boundaries. 
Properties are then distributed across the layer in accordance with the simple 
cubic relations 



where n ~ (y - y^) / (>'2 - y^) • The shear layer is centered about the plume 
interface and grid points are evenly spaced across it in increments of 
A* - (+2 1) . 

4-7.3 B oundary layer o ption .- In the jet entrainment problem, the dis- 
placement effect of the external and nozzle boundary layers must be accurately 
accounted for in the initialization process as will be evident from the 
discussion presented in the next section. In this initialization option, 
schematically depicted in Figure 4,6, the user is required to supply values 
of the displacement thicknesses <5* 2 and nondimensional frictional veloci- 
ties u l,2/ u l,2 at nozzle lip.’ The velocity profile used in the boundary 
layers is given by (see Reference 19, pp, 690-707) 


U 1,2~ U (Q 



2.5 In £ + 1.38 



where Cole's universal wake function 


V (C) 


is well approximated by 


(33) 


w <o = 1 + sin 

Here, £ is the nondimensional variable 

£ = |r - r c |/6 1>2 


(34) 


(35) 


and the frictional velocities u* 2 are related to the wall shear stress by 


U l,2 _ Xw l,2 _ 1 / 9u\ 

u i 2 (p 1 - 1 )^ 2 b^i 2 W/ 


(36) 
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where Re is the Reynolds number per unit length. A nominal value of 
u*/ui 2 c 1/30 is built into the code as a default option for users not sup- 
plying this information. 



program determines 

r i,2 = r C i ^1,2 

using velocity profile 

^4-^ s “2.5AiC~l.38[2w (e) ] 

U |,2 

where 6 = | r — r c | / 2 8 , t2 

Figure 4.6 Boundary layer initialization. 

The temperature variation across the boundary layers is given by the 
Crocco relation, which for a Prandtl number of unity and insulated walls, is 
given by 
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P ■ 1 4 M f , fx . h!1U] 

r l,2 2 1,2 ^ u l,2_ 


If the boundary layer thickness is small in comparison with the nozzle exit 
radius, 6^ 2 are readily obtained by the relations 


iLi = f r _ uoi t i,2 ~] 

S l,2 J n L U l,2 T(S) J 


For thick boundary layers, we start with the expression 




r j +6 l,2 


l.-iS-lLl rdr 

u, - T 


and after some manipulation, obtain 


*1,2 - 1 ! + Vjj + 2I 2 (sjVri + ^.z/ri 




4.7.4 ke 2 initialization .- In the absence of any information regarding 
initial turbulence levels, the initial profiles of k and e are obtained 
from the mixing length model based upon the rel evan t mean flow profiles. The 
ML model relates the turbulent shear stress p u ! v ! to the local velocity 
gradient through the relation 


p u 'v' = p V 


Y9u_\ 3u 

\3y/ 3y 


With the estimate that k = | u f v ’ | /0 . 3 , 


U t (y) |3u/3y[ 
0 . 3p (y) 
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and assuming that P/e = 1 , we obtain 


e( y) = 0.09p(y)k 2 (y). (43) 

v t (y) 

5. CODE VALIDATION 
5.1 Incompressible Free Shear Layers 

The validation process for BOAT was initiated with the analysis of two- 
dimensional incompressible free shear layers. Any code must demonstrate its 
ability to correctly predict this basic flow before being utilized in more 
complex situations. In performing these shear layer calculations, the ratio 
of l/6 in the Prandtl mixing length model was varied and comparisons were 
made with experimental data to "calibrate” the model. The experimental data 
are generally given In terms of the spreading parameter, a , as defined by 
the relation 


1.855 Ax 

< A y).i -.9 


(44) 


where (Ay) i -t- 9 is the change in shear layer width over the distance Ax ; 
the shear" layer ‘width is defined as the distance between the points where 
(u - u 2 )/(u 1 - u 2 ) ■ -1 and .9 . This nomenclature is illustrated in 

Figure 5.1. 



29 






The data were compiled and reviewed by Rodi 20 , the relatively large spread in 
the data being attributed to varying turbulence levels in the free streams. 

The higher values of o (smaller spread rates) correspond to the smaller 
free stream turbulence levels. Since the shear layer boundary conditions are 
based upon zero free stream turbulence, the predictions should match these 
higher values. Figure 5.2 shows that both the ke2 and ML model (with fc/6 
= .065) predictions are in good agreement with these data. In terms of pro- 
file shapes, both the ke2 and ML models adequately predict both the velocity 
and turbulent shear stress profiles as illustrated in Figures 5.3a and b . 

5.2 Incompressible Free Shear Layers with Initial Boundary Layers 

The ability of BOAT to analyze an initially nonsimilar mixing region^has^ 
been assessed by performing Test Case 4 of the NASA Shear Flow Conference ’ • 

In this case, two streams are initially separated by a symmetric airfoil with 
a 10° trailing edge, and initial profile data is supplied at a station 0.1 cm 
downstream of the trailing edge. The dual length scale approach in the ML 
model is required in this situation. The initial velocity profile data are 
depicted in Figure 5.4 (a), and comparisons of BOAT code velocity profile pre- 
dictions with the data at a number of downstream stations are presented in 
Figures 5.4(b) through 5.4(h), respectively. The BOAT calculation was per- 
formed with a 31-point grid and the ratio of 1/6 in the Prandtl mixing length 
model set equal to 0.65. The results in the nonsimilar initial region (i.e., 
x < 25 cm ) are considered to be reasonable while those downstream are quite 
good . 



Figure 5.3a Comparison between predicted (BOAT) and measured velocity 
an d shear stress profile for two-dimensional shear layers, 
ke2 model. 
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Figure 5.3b Comparison between predicted (BOAT) and measured velocity 
and shear stress profile for two-dimensional shear layers, 
Prandtl mixing length model. 

In Figure 5.5, the predicted velocity profiles obtained with BOAT em- 
ploying both the ML and ke2 models are compared with the Lee data at 12.7 cm 
as well as with ML and ke2 predictions by Launder et al. 15 and ML predictions 
by Rudy and Bushnell 22 . The BOAT ML predictions and both ke2 predictions fit 
the data quite well. The ML predictions of both Launder and Rudy do not pick 
up the velocity-defect region at all. Similar trends regarding the predictive 
capability of each model are observed in Figure 5.6, which compares predicted 
and measured shear stress profiles. At x = 76 cm , where the mixing takes 
on a shear layer type of behavior, the BOAT ML predictions (Fig. 5.7a) are 
quite good in contrast to the predictions of the ML models in the other two 
codes. The BOAT kc2 predictions (Fig, 5.7b) are also quite good while those 
in the GENMIX code are somewhat poorer. 

These comparisons show that both the ML and ke2 models in BOAT accurately 
predict a type of flow situation, which is quite analogous to the nearfield 
mixing layer downstream of nozzle boattails. In addition, they demonstrate 
that the performance of a turbulence model cannot be judged separately from 
the code in which it is contained. The noted differences in the model/code 
predictions are attributable to both minor variations in the specific model 
formulations within the codes as well as to variations in the computational 
aspects of the various codes such as grid distribution, shear layer growth 
expression, etc. 
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Figure 5.6 Comparison between predicted and measured shear stress profiles 
for 2D shear layer with initial boundary layers. 
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Figure 5.7b Comparison between predicted and measured velocity profiles 
for 2D shear layer with initial boundary layers. 

Although the above comparisons are for incompressible flows, we expect 
that the BOAT code will be equally valid for the mildly compressible flow 
situations in typical aircraft exhausts (i.e., jet exhaust Mach numbers of 
approximately one, subsonic/transonic external streams). As demonstrated in 
the compilation of jet spreading rate data presented in Reference 13, notice- 
able reductions due to compressibility are associated with jet Mach numbers 
greater than 1.5 to 2 . 
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5.3 Coaxial Jet Mixing 


Although nearfield mixing is of prime concern in this study, it is also 
of interest to demonstrate the validity of BOAT for farfield situations. This 
has been done via comparisons with coaxial jet mixing data. The value of 
£/6 = .08 used for farfield axisymmetr ic situations in the ML model was, in 
fact, arrived at by matching farfield decay rates for jets into still air. In 
application of the ML model, the value of the ratio of £/6 is changed from 
.065 to .08 abruptly when the mixing layer reaches the axis. It is expected 
that reasonable rates of mixing will be predicted for both the nearfield shear 
layer and in the farfield fully developed mixing layer. However, poorer 
agreement is anticipated in the transitional region where the turbulence 
processes are well out of equil ibr ium , and simple models predicated upon equi 
librium assumptions may not be adequate. These statements are confirmed by 
the predictions achieved with the ML model in BOAT in analyzing the Forstall 
and Shapiro data 23 (Test Case 9 of the NASA Shear Flow Conference) for coaxial 
jet mixing. Figure 5.8 indicates that while the rate of velocity decay in the 
farfield is well predicted by the ML model, the performance in the transi- 
tional region (5 < X/Dj < 30) is poor and affects the overall quality of the 
results in the farfield. The ke2 model, however, does admirably in analyzing 
the nearfield, farfield, and transitional regions. 

AXISYMMETRIC JET IN MOVING STREAM 



Figure 5.8 Comparison between centerline velocity decay 
predictions and data for an axisymmetric jet. 
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6. APPLICATIONS TO LABORATORY COLD AIR JETS 


6.1 Displacement Thickness Type Correction to Inviscid Plume Geometry 

For the cases to be considered below (which contain many of the essential 
features of actual aircraft exhaust flows), the external boundary layer thick- 
nesses have been rather substantial (i.e., comparable in size to the nozzle 
exit radius as depicted in Figure 6.1). Thus, the streamline deflections 
associated with the entrainment process are primarily induced by the mixing 
away of the low velocity (wake-like) initial region. These deflections are 



Figure 6.1 Determination of effective inviscid plume boundary 
shape for large initial boundary layers. 

comparable In magnitude to (or even larger than) those produced by purely 
inviscid processes. For analyzing the flow In such situations, the following 
procedure has been followed. In view of the rather small deflection angles 
along the plume interface, mixing processes along the interface are analyzed 
In standard cylindrical coordinates rather than in plume-oriented boundary 
layer coordinates. Then, the normal velocity, v e (x) , at the outer shear 
layer edge is in the radial direction and includes contributions from both 
viscous (entrainment) and inviscid (blockage) processes. To account for the 
radial variation in v required to satisfy the continuity equation, the 
Injection velocity, vj (x) , to be applied along the interface for the sub- 
sequent external flowfield calculation, is obtained by modifying v e (x) via 
the relation 


v j (x) 


P e (x)r e (x) 

" p J r J 


V e (x) 


(45) 
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Then, the "effective” plume boundary, r eff (x) , is determined by adding 
6*(x) (obtained from d6*/dx = vj/u e ) to r e Tj , the 6*(x) integration 
being initialized with the value of 6* at the nozzle exit from the boundary 
layer calculation. 

In terms of incorporating BOAT into the iterative patched LRC method- 
ology, the following sequence of calculations is entailed: 

(1) Iterative solution of exhaust/external flows without entrainment to 
yield inviscid flowfield maps and boundary layer properties at 
nozzle exit. 

(2) Overlaid shear layer calculation employing inviscid flowfield maps 
to yield streamline deflections due to entrainment process. 

(3) Plume geometry modification accounting for streamline deflections 
due to entrainment* 

(4) External flow calculation repeated over afterbody and modified plume 
geometry. 

(5) Exhaust flow calculation repeated employing pressure variation 
determined from external flow calculation. 

Steps (2) through (5) can be repeated in a successive fashion until a con- 
verged viscous/inviscid flowfield solution is obtained. In this approach, 
the entrainment effect directly modifies the external flow structure and 
indirectly modifies the exhaust plume structure. The optimal procedure to 
be utilized in arriving at an overall converged solution has not been formal- 
ized. An assessment of alternative procedures is currently being investi- 
gated by R.G. Wilmoth at LRC. The cases reported below do not represent 
fully converged solutions but rather first-order entrainment corrected solu- 
tions employing steps (1) - (4) above. 

6.2 Effective Cold Air Jet Boundaries 

Preliminary results using BOAT within the patched LRC methodology have 
been obtained for one of the boattail nozzle configurations tested with cold 
air jets by Reubush 12 . Experimental boattail surface pressures are presented 
in Reference 12 for various free stream Mach numbers and jet total pressure 
ratios (Pt j/P e ) • BOAT calculations were performed for a fully expanded 
(P t =,!? e = 2) and a modestly underexpanded jet (Pt,j/ p e “ 3) with a nominal 
velocity ratio, Wu, , of about 0.5 , and free stream Mach number of 0.4 . 
Initial external boundary layers at the end of the boattail were quite thick 
in both cases. 

The effective boundaries downstream of the nozzle exit for the two cases 
investigated are shown in Figure 6.2 ( x/r. - 0 corresponds to the end of the 
boattail, which is also the nozzle exit). Also shown are the inviscid plume 
interfaces obtained by the method of Reference 5 and the effective solid 
plume boundaries obtained by treating the inviscid plume interface as a solid 
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body to which is added the calculated boundary layer displacemnt thickness. 
This latter boundary is commonly used to predict jet plume blockage effects 
without considering entrainment (e.g., see Reference 2). The important 
observation here is that when the effects of entrainment are considered, the 
effective plume boundary is substantially reduced in size primarily because 
the mass deficit, associated with the initial boundary layers, mixes away 
rather rapidly. Thus, treating the inviscid plume as a solid boundary is 

INVISCID PLUME 

EFFECTIVE SOLID PLUME 

(BOUNDARY LAYER ADDED) 



Figure 6.2 Comparison between effective boundary calculated by BOAT and 
effective solid plume for fully expanded and underexpanded 
jet mixing, M e * 0.40 , ke2 turbulence model. 
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quite unrealistic. Since jet aircraft nozzles are generally located in 
regions of thick boundary layers (e.g., near the aft end of the aircraft), a 
similar reduction in effective boundary size is expected, 

6.3 Predicted and Measured Boattail Pressures 

Boattail pressure distributions predicted by the inviscid code of Refer- 
ence 4 using the effective solid plume (no entrainment) and the effective 
boundary (with entrainment) are compared with the experimental data in 
Figures 6.3a and b. A sketch of the boattail geometry and inviscid plume 
shape is also indicated. Excellent agreement with the experimental pressure 
distribution is obtained at both jet pressure ratios when entrainment effects 
are included, with a corresponding improvement in the predicted boattail drag 
coefficient, C D g . While the effect of underexpansion on the inviscid 
plume geometry appears relatively small (see Fig. 6.2), the pressure distri- 
bution is quite sensitive to this effect as evidenced by the variations in 
Cp downstream of the nozzle exit (Fig. 6.3). Including the (viscous) effect 
of entrainment decreases these variations, indicating a weakening of the 
effective wavelike amplitude of the inviscid plume. Such an effect is quali- 
tatively similar to that produced by the boundary layer on the flow about a 



Figure 6.3a Effects of jet entrainment on boattail pressure 
distributions and comparison with experimental 
data, Me = 0.40 , ke2 turbulence model. 
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compression corner or expansion shoulder, i.e., a weakening of the "inviscid" 
wave structure. Based on these limited comparisons, the interactive entrain- 
ment model appears to provide a reasonable representation of the boattail 
f lowf ield . 



Figure 6.3b Effects of jet entrainment on boattail pressure 
distributions and comparison with experimental 
data, M e = 0.40 , ke2 turbulence model. 

6. A Sensitivity Studies 

By varying some of the input and modeling parameters required for en- 
trainment calculations, the relative dependence of the effective geometry 
(and hence boattail pressures) on these parameters can be established. Areas 
explored were initial boundary layer properties, turbulence models, and 
pressure gradients. 

6.4.1 Sensitivity to initial boundary l ay er properties .- In the cases 
studied here (and for most cases of practical interest), the internal (nozzle) 
boundary layers are small with respect to the nozzle exit radius while the 
external (afterbody) boundary layers are relatively large. For the boattail 
configuration used for the calculations, nominal values of 6*/r j - .013 
(taken from Yaros 10 ) and 6*/*^ *24 were employed. Sensitivity studies 

in which the size (6*) and shape (u*) parameters of the internal boundary 
layer were varied by a factor of 2 about their nominal values produced 
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negligible changes in the effective plume geometry, as anticipated in view of 
the negligible internal boundary layer size. A similar variation of the shape 
factor for the external boundary layer (nominal value of U 2 ~ .033), holding 
6 ? fixed, also produced negligible changes. Substantial variations in 62 
will, of course, have an appreciable effect on the effective plume geometry. 
However, since 6 * is supplied as part of the iterative external flow solu- 
tion, it cannot be arbitrarily varied without further iterations between BOAT 
and the external flow code. This is quite time consuming with the present 
manual coupling between codes and not warranted in view of the rather obvious 
sensitivity of the effective geometry upon 6 * . 

6.4.2 Sensitivi ty to turbulence models .- By performing the calculations 
with both a simple (ML) and more detailed (ke2) turbulence model, the level 
of sophistication required in the turbulence modeling can be assessed. Com- 
parisons of the predicted effective boundary shape for both cases studied are 
given in Figure 6.4. In both cases, the ML model predicts a more rapid rate 



Figure 6.4 Effect of turbulence model on effective plume boundary shape. 
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of mixing (and corresponding rate of entrainment) than the ke2 model and thus 
a narrower effective plume geometry. The difference is substantially more 
pronounced in the underexpanded case where appreciable pressure gradients 
exist and the shear layer edge conditions are thus varying. Inviscid calcu- 
lations over these geometries (and those for other similar test cases) have 
indicated that either the ML or ke2 model can yield a reasonable estimate of 
entrainment effects for fully expanded jets while the ke2 model produces more 
consistent results for underexpanded jets. 

A. 3 Sensitivity to pressure gradients ,- To test the sensitivity of 
the present results to pressure gradients, BOAT calculations (with the kc2 
model) were repeated with the pressure gradient terms in both the momentum 
and energy equations set identically to zero (while local values of density 
were computed employing the mapped inviscid pressures). A comparison of the 
effective plume geometry with and without pressure gradients is shown in 
Figure 6.5. With the pressure gradients deleted, the entrainment effect is 
markedly reduced; in fact, the effective boundary shapes are only slightly 
improved from those obtained treating the plume interface as a solid boundary 
(see Fig. 6.2 for comparison). Clearly, the favorable nearfield pressure 
gradients contribute strongly to accelerating the flow in the low velocity 
mass-deficit region and thus increase the overall entrainment rate. This 
marked effect of pressure gradients clearly rules out the isobaric assump- 
tions used in previous entrainemnt models 10 ’ 11 . 

7 . CONCLUDING REMARKS 

1. The use of a "displacement thickness" (or normal injection velocity) 
correction to the inviscid plume boundary has been introduced to account for 
the effects of jet entrainment on the inviscid external subsonic/ transonic 
flow calculation. A computational model (BOAT) has been developed which pre- 
dicts the rate of jet entrainment via an overlaid, parabolic procedure from 
which the displacement thickness can be determined. Limited comparisons 
between predicted and measured boattail pressure distributions have been 
quite favorable, indicating that use of the BOAT code to calculate jet 
entrainment, in conjunction with the patched NASA/LRC system, shows great 
promise as a computational approach for predicting nozzle boattail drag. 

2. The BOAT code has been shown to be a viable tool for calculating 
nearfield turbulent mixing processes. The application of both eddy viscosity 
(Prandtl mixing length) and two-equation (ke2) turbulence models within BOAT, 
enable it to adequately predict basic (constant pressure, nonreacting) free 
shear layers. Sensitivity studies have indicated that the ke2 model, which 
accounts for the turbulence "history," yields overall results of better 
quality for underexpanded, variable pressure plumes. The assessment of 
BOAT's capabilities at higher Mach numbers and for reacting/high temperature 
exhausts via comparisons with a broad based body of data (supersonic jets/ 
shear layers, turbulent diffusion flames, rocket exhaust plume data, etc.) is 
in progress. 

3. The boattail pressure changes produced by the entrainment correction 
to the inviscid plume boundary suggest that the weak viscous/inviscid 
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interaction approach adopted for this study may also be adequate for more 
generalized situations. The rather large influence of pressure gradients on 
the entrainment rate, in application of the overlaid methodology, clearly 
demonstrates the inadequacy of isobaric mixing assumptions employed in 
earlier modeling efforts. Further studies are in progress to assess the 
adequacy of the viscous-inviscid coupling procedure over a wider range of 
operating conditions. 

4. An automated procedure for coupling the results of the various codes 
comprising the patched viscous/inviscid model is in progress at LRC . Toward 
this end, a new inviscid exhaust plume code has been developed which will be 
directly coupled with the BOAT code providing a one-pass solution of the 
plume and mixing layer. The inviscid code employs shock— capturing method- 
ology 24 ’ 25 in a mapped coordinate grid and has been tailored for ease of 
operation and extremely fast run time. 



Figure 6.5 Effect of pressure gradient on effective plume 
boundary shape, ke2 turbulence model. 
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5« An evaluation of turbulence models for application to afterburning 
rocket and aircraft plumes is in progress (under NWC support) using the BOAT 
code. Specific objectives and preliminary findings are reported in Reference 
26. The study should yield valuable information regarding the choice of tur- 
bulence models in extending BOAT to situations with combustion and marked 
compressibility. 
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